function hasil = heun(f,x0,y0,b,h)
%HEUN fungsi menghitung nilai turunan dengan metode Heun
%   
    n = (b-x0) /h;
    y = y0;
    x = x0;
    
    syms a b    
    for r=1:n
        ys = y;
        y = y + h*subs(f,[a b],[x y]);                       %PREDICTOR dengan Euler
        y = ys + h/2 * (subs(f,[a b],[x ys]) + subs(f,[a b],[x+h y])); %CORRECTOR dengan Heun
        
        x = x + h;
    end
    hasil = y;   
end

